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A study of donor electron spins and spin-dependent electronic transitions involving phosphorous ( 31 P) atoms 
in proximity of the (111) oriented crystalline silicon (c-Si) to silicon dioxide (SiCh) interface is presented for 
[ 31 P] = 10 15 cm" 3 and [ 31 P] = 10 16 cm" 3 at about liquid 4 He temperatures (T = 5 K- 15 K). Using pulsed 
electrically detected magnetic resonance (pEDMR), spin-dependent transitions between the 31 P donor state 
and two distinguishable interface states are observed, namely (i) P(, centers which can be identified by their 
characteristic anisotropy and (ii) a more isotropic center which is attributed to E' defects of the SiC>2 bulk 
close to the interface. Correlation measurements of the dynamics of spin-dependent recombination confirm 
that previously proposed transitions between 31 P and the interface defects take place. The influence of these 
electronic near-interface transitions on the 31 P donor spin coherence time T% as well as the donor spin-lattice 
relaxation time 7j is then investigated by comparison of spin Hahn-echo decay measurements obtained from 
conventional bulk sensitive pulsed electron paramagnetic resonance and surface sensitive pEDMR, as well as 
surface sensitive electrically detected inversion recovery experiments. The measurements reveal that both Ti 
and T\ of 31 P donor electrons spins in proximity of energetically lower interface states at T < 13 K are reduced 
by several orders of magnitude. 
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I. INTRODUCTION 

Due to their technological importance, the properties of 
phosphorous ( 31 P) donors in crystalline silicon have been in- 
vestigated extensively for more than half a century. During 
this time, magnetic resonance based methods have revealed 
many aspects of the microscopic nature of 31 P as well as the 
electronic processes in which it is involved. Since the first 
mapping of a donor wavefunction using an electron nuclear 
double resonance technique 1 , a large and diverse amount of 
information regarding the electronic and quantum mechanical 
properties of such systems has been reported 2 . This wealth of 
information has contributed to the technological exploitation 
of silicon to the extent that it has become the most widely uti- 
lized semiconductor in the global electronics industry. How- 
ever, even with all the information regarding silicon available 
to us, there are still questions regarding the ability to exploit 
the quantum mechanical nature (specifically spin) of dopants 
and charge carriers for technological applications such as spin 
transistors 3 or quantum computers 4 , the latter of which are 
proposed to utilize electron or nuclear spins of phosphorus 
donors in silicon as quantum bits. While there has been sig- 
nificant experimental effort for the implementation of these 
and other concepts 5 , challenges remain in a number of differ- 
ent areas, the most pressing being the difficulty in detecting 
the spin of individual donors without reducing its quantum 
memory time (the coherence time, 72) 6 . Nevertheless, donors 
in silicon retain promise in this area due to their extremely 
long phase coherence times, with T2 exceeding 60 ms for the 
electron spin 7 - and 1 second for the nuclear spin 8 at liquid He 
temperatures. 

Many of the recent schemes for silicon based spin- or quan- 



tum devices involve electronic processes occurring at or near 
interfaces, particularly the c-Si/SiC>2 interface. This presents 
advantages since locating spins near interfaces allows them to 
be controlled with surface gates 4 , and detected with surface 
electronics 9 ' 1 °i 1 1 ' 1 2 ' 1 3 . However, it may also lead to a decrease 
in spin coherence due to the spin-spin interactions with sur- 
face state s 14 ' 15 , as well as to the loss of quantum informa- 
tion following spin-dependent recombination through surface 
states. Since it is not possible to obtain defect free c-Si/SiC>2 
interfaces^ (we note that work is currently being undertaken 
to obtain defect free interfaces using H-termination but it is in 
the very early stages 1718 ) and since there are even proposals 
to use interface defects as probe spins to readout 9 the spin of a 
single donor, the understanding of interface electron and spin 
transitions has become important. 

In the following, pulsed electrically detected magnetic reso- 
nance (pEDMR) and pulsed electron paramagnetic resonance 
(pEPR) measurements performed on (111) surface oriented 
c-Si samples with [P]=10 15 cm" 3 and [P]=10 16 cm" 3 at tem- 
peratures between T = 5 K and T = 15 K are presented. The 
data was collected in order to elucidate the nature of spin- 
dependent electronic interface transitions involving 31 P donor 
and interface states and to then determine how these processes 
influence the coherence time of the 31 P donor electron spins 
in proximity of these interface defects. The latter was ac- 
complished by comparison of the interface sensitive pEDMR 
measurements to bulk sensitive pEPR measurements of 31 P 
donors. Our results are discussed with regard to their implica- 
tions for the ability of spin qubit readout using interface defect 
probe spins whilst maintaining long coherence times. We em- 
phasized that while the key questions motivating this study are 
centered about 31 P qubit coherence times, the study presented 
follows an extensive number of previous EPR 1 9,2Q ' 2 1 i 22 - 23 i 24 i 25 ; 
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EDMR 10 ' 25 ' 26 - 27 ! 28 and pEDMR 29 i 3Q i 31 studies carried out on 
various c-Si/Si02 interface defects as well as electronic trap- 
ping and recombination processes of interfaces with different 
surface orientations and 31 P doping concentrations. Most of 
these studies aimed to enhance the understanding of electronic 
processes relevant for materials systems used in conventional 
c-Si based microelectronics and photovoltaic devices. Thus, 
the study presented here may also be of relevance for conven- 
tional silicon technologies. 



II. EXPERIMENTS 

Following the recent demonstration of electrical detection 
of 31 P spin states 10 we anticipate, for the pEDMR measure- 
ments presented here, the presence of spin-dependent recom- 
bination between 31 P donor states and energetically lower c- 
Si/SiC>2 interface states, as illustrated for the example of the 
Pb interface state in Fig.Q] FigureQla) displays a geometrical 
stick and ball sketch of the c-Si/SiC>2 interface where some 
of the silicon (111) surface atoms posses unsaturated bonds 
(called dangling bonds). These dangling bonds are all ori- 
ented along the (111) axis 2 ^ and are highly localized states 
with much p-content and only small s-conten t 32 ' 33 . Silicon 
dangling bonds oriented along the (111) direction (also called 
Pb centers) are highly anisotropic and, since they are para- 
magnetic when uncharge d 19 ! 34 , can be identified using EPR 
by their well known Lande (g) factor anisotropy 25 , revealed 
by measurement of EPR spectra as a function of the angle 
between the Z?o field and the (111) direction (see sketch). 
In contrast to the Pb state, the shallow 31 P donor state (illus- 
trated by the large shaded circle in Fig. [TJ a )) is a significantly 
less localized state having an s-like envelope function with a 
localization length of » 3nm, strongly modulated by the pe- 
riodic crystalline structure of its host environment 35,36 . Since 
the 31 P donor state encompasses thousands of silicon sites, an 
exchange interaction between near-surface 31 P donors and Pb 
centers can occur allowing the transition of the donor electron 
into the energetically lower interface state. This transition is 
sketched for an interface band diagram in Fig. Q2b)- The il- 
lustration shows that: (i) Both the 31 P donor state as well as 
the Pb state become charged (positively and negatively, re- 
spectively) through the transition, (ii) Because of the Pauli 
principle, and spin conservation due to the weak spin-orbital 
coupling of silicon, the transition is spin-dependent. Hence, 
the transition provides a spin-to-charge conversion mecha- 
nism and as a result, since the charge of the 31 P + /Pb ~ can be 
detected through the measurement of recombination currents, 
jl P /Pb pairs have been proposed as electric readout mecha- 
nism for 31 P qubits 9 *i£. 

Experimentally, time domain measurements of spin- 
dependent c-Si/SiC>2 -interface recombination were con- 
ducted in order to verify that the qualitative behavior of in- 
terface currents after EPR excitation of the 31 P or the Pb states 
is similar to higher doped c-Si samples with (100) orientation. 
In order to also verify electrically detect spin-Rabi oscillation 
(as expected from previous studies&i^), transient-nutation- 
like measurements of the integrated current transients as a 
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Figure 1: (color online) (a) Ball and stick illustration of the atom 
scale structure of the c-Si/Si02 interface. The 31 P donor state (green 
circle) is significantly larger than the Pb center. Exchange coupling 
between the two states is possible when they are sufficiently close to 
each other. The orientation of the interface is defined by the angle 9 
between the externally applied Bo magnetic field and the (111) crys- 
tal axis which is perpendicular to the crystal surface, (b) Band sketch 
illustrating the spin-dependent 31 P /Pj, charge carrier recombination 
transition. For details see text 



function of the excitation lengths were carried out. Follow- 
ing these preparatory measurements, four experiments were 
performed to address the following questions: 

(i) Is the EDMR signal that has previously been observed 10 
at magnetic fields between the two hyperfine split lines of the 
31 P for measurements conducted at X-band excitations (ap- 
prox. 10 GHz) truly due to Pb states? If so, are Pb states 
the only interface states involved in spin-dependent transi- 
tions as illustrated in Fig. fflb) or are there other states at or 
near the interface, or even in the bulk, which could contribute 
to the observed signals? In order to address this question, a 
systematic study of the EDMR spectrum (the magnetic field 
dependence) as a function of the interface orientation angle 9 
was made in order to observe all ^-factors involved in spin- 
dependent recombination and to detect possible anisotropies 
of these centers. 

(ii) Are the observed spin-dependent processes truly due to 
pairs of 31 P and interface states as depicted in Fig.[TJb) or are 
the signals corresponding to different ^-factors due to inde- 
pendent processes? In order to address this question, a series 
of different spin-dependent current transients was recorded 
under various temperatures, light intensities, surface orienta- 
tions, sample voltages and offset currents as well as excitation 
powers and lengths. For all applied conditions, the current 
transients were recorded after (a) resonant excitation of the 
two 31 P lines and (b) the interface states and the correlation 
between the dynamics of 31 P and interface states were com- 
pared. 

(iii) As an understanding about the nature of the interface 
defects and the transitions between interface defects and the 
31 P donor electrons was established, the main question of this 
study was addressed: Are the coherence times of 31 P donors 
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near interface defects compromised? For this, a series of co- 
herence time (T2) measurements was carried out on 31 P donors 
in proximity to interface states and with the same samples, 
under the same conditions (in fact, during the same experi- 
mental runs), the coherence time of bulk 31 P was measured 
for comparison. For the interface T2 measurements, modified 
Hahn echoes, detected with pEDMR, were used, in a similar 
way as recent studies of 31 P doped c-Si samples with (100) 
orientation 6 and spin-dependent 31 P bulk processes at very 
high magnetic fields 37 . For the bulk T2 measurements, con- 
ventional ESR detected Hahn-echo experiments were carried 
out. This comparative study of T2 times was made as a func- 
tion of the temperature for 5 K < T < 13 K. 

(iv) Finally, in order to obtain information on wether the 
spin-coherence times T2 of near-interface defect 31 P is deter- 
mined by the electronic transitions between 31 P and interface 
states or by the interface defect induced spin-spin relaxation 
processes as suggest by de Sousa 15 , a comparison of elec- 
trically detected T2 times and T\ times was made for a tem- 
perature of T = 5L. In order to measure T\ , electrically de- 
tected inversion recovery measurements were performed. To 
the knowledge of the authors, this is the first time that such an 
experiment is demonstrated. 



HI. EXPERIMENTAL DETAILS 

For the experiments presented in the following, we used 
300 /xm thick, (111) surface oriented silicon since Pb centers 
for this surface orientation are all identically oriented. This 
is in contrast to the (100) surface previously used for similar 
experiments 610 as the (100) surface allows silicon dangling 
bonds (here they are called Pbo states) to exist in two orienta- 
tions at the same time which makes their mutual EPR spectro- 
scopic distinction as well as their distinction from other pos- 
sible interface defects very difficult. The experiments were 
conducted with dopant concentrations of [P] = 10 15 cm~ 3 and 
10 16 cm~ 3 which is less than 10 17 cm~ 3 of previous EDMR 
studies 610 . The lower concentrations ensured that interactions 
between neighboring 31 P atoms which are known to exist 3 ^ at 
[P] = 10 17 cm~ 3 can be neglected. 

In order to enable the electrical detection of spin- 
coherence, the c-Si samples required contact structures which 
would not strongly distort the homogeneity of the spin res- 
onant microwave fields B\. This problem was solved using 
a method similar to previous pEDMR studies at X-ban d 10 ' 39 
- we designed long match-like sample substrates on which 
electrical sample contacts outside of the microwave field are 
connected to the sample via ~ 50 mm long thin-film Al wires 
whose thickness of 100 nm is below the penetration depth of 
the applied microwave radiation. The contact structures were 
fabricated using a photo lithographical lift-off process that 
was carried out after the Al film was evaporated on the H- 
terminated silicon (111) surface of the c-Si samples that had 
been prepared by a wet treatment with hydrofluoric acid. In 
order to maximize the density of interface states, a native ox- 
ide was grown on the (111) surface after the contact deposition 
by exposure of the sample to ambient air. 



All EDMR and EPR experiments were carried out at X- 
band using a cylindrical dielectric low-Q pulse resonator 
which was part of a Bruker Elexsys E580 EPR spectrometer. 
The sample temperatures were obtained with a 4 He flow cryo- 
stat, the excess charge carriers were induced through a spectral 
cold light source (IR and UV filtered spectral light) with an 
incandescent light source (Schott KL 2500 LCD) with lamp 
temperature of about 3000 K producing an integrated spec- 
tral intensity of approximately 5 Wcm~ 2 at the sample sur- 
face. The EDMR experiments were conducted by establishing 
a constant offset photocurrent using a constant current source 
with a time constant in excess of the experimental shot repeti- 
tion time. Current transients were digitized and recorded fol- 
lowing current amplification by a Stanford Research SR570. 

The raw data recorded for the presented measurements was 
a combination of spin-dependent currents and microwave in- 
duced artifact currents. The latter can be recorded separately 
by measurement of the current response at off-spin resonant 
Bo-fields. Magnetoresistance effects on the microwave in- 
duced currents can be linearly extrapolated for c-Si at the 
given magnetic fields. The microwave current transients ob- 
tained from this procedure were subtracted from the raw data 
in order to reveal the current transients solely caused by spin- 
dependent transitions. 

IV. MEASUREMENT RESULTS 
A. Identification of spin-dependent transitions 

1. Experimental data 

In order to confirm the results of Stegner et al^ for a c- 
Si:P/SiO? interface with smaller 31 P concentration and a (1 1 1) 
surface orientation of the silicon sample, transient measure- 
ments of photocurrent changes Al were recorded under vari- 
ous illumination conditions and temperatures. Figure [2] dis- 
plays a data set of Al (Bo,t) recorded as a function of the mag- 
netic field Bq and the time t after a 96 ns long microwave pulse 
with a frequency of / = 9.749 GHz and a power of P sa 8 W 
for T = 5 K and a constant photocurrent of / = 270 jJ.A. The 
sample orientation was 0=0. The data set clearly confirms 
the expected EPR induced currents with three local response 
maxima at B Q = 346.37 mT, B = 347.9 mT and B = 350.55 
mT. The two outer peaks which are separated by a magnetic 
field of « 4.2 mT are the two hyperfine lines of the 31 P donor 
electron while the peak close to the low field 31 P line has been 
attributed to interface defect states 10 . 



2. Discussion 

The transient behavior at the magnetic fields with EPR re- 
sponses confirms the measurements by Stegner et alJ^. It 
consists of a brief photocurrent quenching after the pulse, at- 
tributed to an enhancement of the interface recombination, 
followed by a longer-lived current enhancement. This en- 
hancement arises because the singlet recombination rate re- 
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turns to its steady state faster than the resonantly quenched 
triplet recombination rate, causing a net quenching of the re- 
combination rat o 10 ' 40 . Note that variations of the tempera- 
ture, the sample voltage and therefore the photocurrent as 
well as the illumination conditions change the quantitative 
dynamics of the observed transients - they do not, however, 
change the qualitative quenching/ enhancement behavior dis- 
played by the data set of Fig. [2] which were also observed 
in previously reported measurements^ on c-Si(100):P/SiC>2 
with [P] = 10 n cm- 3 . 

Both the quenching and the enhancement shown in Fig.[2|d) 
are well fit with simple exponential decay functions, in agree- 
ment with the observation of Stegner et al 10 . This is somewhat 
counterintuitive since the random spatial distribution of the 
paramagnetic states involved in these transitions (see discus- 
sion in section HV Bl ) suggests that the distances between pairs 
of paramagnetic states between which spin-dependent transi- 
tions occur are widely distributed. This distance distribution 
of states also implies a distribution of transition times 41 which 
suggests that the observed current transients should be distri- 
butions of exponential functions for which fits with single ex- 
ponentials would be poor. In contrast, observed transients are 
in good agreement with single exponential fits which suggests 
that only a narrow range of transition times exist and there- 
fore, only pairs with a narrow range of intra-pair distances 
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Figure 2: (color online) Plots of the microwave pulse (at / = 
0) induced change AI of an interface photocurrent along a c- 
Si(l 1 l):P/SiC>2 interface as a function of time t and the applied mag- 
netic field Bq. (a,c) Plots of AI as functions of So for the times 
/ = 6.1 /is and t = 38 fls, respectively, for which a quenching and 
an enhancement of the current changes reach their extrema. (b) color 
contour plot of the entire data set AI(t,Bo) containing the data of the 
plots in (a), (c) and (d). (d) Plot of AI as a function of the time t for a 
magnetic field So = 350.5 mT which was on resonance with the high 
field 31 PEPR line. 



contribute to the observed signals. We attribute the existence 
of a "main-pair distance" which dominates the observed sig- 
nals to two factors: first, pEDMR signals vanish for pairs with 
very large distances (larger than the localization length of the 
two paramagnetic states 41 ) as the probability for recombiation 
is greatly diminished. Second, for very short distances, where 
the exchange between the two states exceeds the Larmor fre- 
quency difference within the pair 42 , the signal vanishes as the 
resonantly induced change of spin pair states between singlet 
and triplet configurations becomes increasingly forbidden 42 . 
Thus, there will always be a finite main-pair distance where 
spin-resonance induced rate changes become maximal and 
pairs around this distance will dominate the observed spin- 
dependent currents. 



B. Identification of interface defects 

1. Experimental data 

The data presented in Fig. [2] confirms that spin-dependent 
interface recombination processes can be observed with 
pEDMR at the c-Si(l 1 l):P/Si0 2 interface which are qualita- 
tively similar to those seen previously at the c-Si(100):P/SiO2 
interface with higher 31 P concentration. It is therefore pos- 
sible to systematically study the nature and the origin of the 
resonances found. Figure [3fa) displays the magnetic field de- 
pendence of the pulse induced photocurrent change AI (Bq) at 
times t after the pulse when the photocurrent changes were 
maximal (in the following referred to as pEDMR spectra) for 
five different surface orientation angles (0° < 9 < 90") and 
two 31 P concentrations. Note that while AI < 0, the data 
has positive signs as each displayed spectrum was normal- 
ized to its respective extremum. The normalization was made 
for better comparison of the data sets since the signals ob- 
tained from samples with [P] = 10 16 cm~ 3 were significantly 
stronger compared to signals from samples with [P] = 10 15 
cm~ 3 (note the higher relative noise in the latter spectra). 

The fits of the pEDMR spectra required at least four dif- 
ferent Gaussian peaks. A comparison with the 31 P EPR and 
EDMR spectra found in the literature 10 - 43 allows immediate 
identification of the peaks at the highest and lowest magnetic 
fields as the well known hyperfine split 31 P donor electron 
resonances. There are at least two additional resonance sig- 
nals present. As these additional peaks are very close to each 
other and to the low field 31 P hyperfine resonance, significant 
ambiguity for the 12 fit parameters (peak centers, width, mag- 
nitude for all four lines) is present. This was overcome by 
a stepwise fit of the spectra: we first determine the 31 P spec- 
trum and subsequently fit the two remaining non-phosphorous 
lines (eliminating six fit parameters). The separation of the 
low-field 31 P -hyperfine line from the strongly overlapping 
non-phosphorous lines was achieved by first fitting the high 
field 31 P hyperfine peak, which has little or no overlap with 
the other resonances, with a single Gaussian line. From the 
result of this fit, we can determine both the position (assuming 
a hyperfine splitting of A =4.18 mT as verified by bulk EPR 
measurements) and shape (assuming only negligible nuclear 
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polarization, which is justified for the given sample tempera- 
tures and magnetic fields Bq) of the low-field 31 P -hyperfine 
peak. The residue following these fits is then able to be fit 
with two Gaussian resonances. The full fit, as well as the 4 
constituent peaks, is shown for all spectra in Figure[3] 

Figure [3jb) shows the g-factors of the non- 31 P lines as a 
function of the angle 9, obtained from the fit. It also displays 
a solid line which represents literature values for EPRs£ and 
EDMR^ 3 ! detected P b centers. 
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Figure 3: (a) Plots of AI as functions of Bq at arbitrary times t the 
after a microwave pulse with arbitrary length T, frequency / 9.5 
GHz and a power P = 8W and under otherwise identical conditions 
as for the data in Fig. [2] The data was collected for five sample 
orientations 8 and two 31 P concentrations. The displayed plots are 
normalized to the maximum of AI. The solid lines represent fits of 
the data consisting of four Gaussian peaks, two related to the 31 P 
hyperfine line and two peaks related to interface defects. The plot for 
[ 31 P ]= 10 16 cm~ 3 and 8 = 90° has the peak assignments to the low 
field (LF) and high field (HF) 31 P resonances, the P b peak as well as 
peak 2. (b) Plots of the fit results of the g-factors for the two interface 
defect lines as a function of 8 for the two 31 P concentrations. The 
solid lines indicate the literature values for the P b center. 



2. Discussion 

The anisotropy and absolute value of the g-factor of one of 
the two peaks is in excellent agreement with the P b literature 
values for both 31 P concentrations. Based on this agreement, 
this peak can be assigned to spin-dependent transitions (re- 
combination) which involve P b centers. Note that previous 
pEDMR measurements on c-Si(l 1 l)/SiC>2 interfaces with no 
31 P doping have shown P b signal a 29 ' 31 . Since there are no 31 P 
pEDMR lines in intrinsic c-Si, it is clear that spin-dependent 
P b interface recombination does not necessarily require the 
presence of 31 P atoms. Thus, the P b involvement revealed 
by the data in Fig. [3] may either involve 31 P as illustrated in 
Fig. [Tfb) or be due to an independent interface recombination 
process. 

In contrast to the P b resonance line, the assignment of the 
second non- 31 P peak (referred to as peak 2, see peak assign- 
ment in Fig. [3j is less straight forward. The fit results for 
[P] = 10 15 cm~ 3 suggests that peak two is isotropic, or only 
weakly anisotropic, whereas the data for [P] = 10 16 cm~ 3 is 
consistent with an anisotropic peak. There are two explana- 
tions for this difference between the two sets of spectra, (i) that 
the nature of peak 2 is different at higher 31 P concentrations, 
suggesting that the observed processes may be different and 
(ii) that the fit error of the center g-factor of peak 2 becomes 
increasingly inaccurate for larger values of 9. The latter may 
be due to the width of peak 2 and the relative weakness of the 
peak 2 intensity in comparison to the signal strengths of the P b 
and 31 P EDMR signals at higher dopant concentrations, be- 
coming increasingly problematic for the fit accuracy at large 
9 as three resonances (low field 31 P , P b , peak 2) have larger 
overlap at higher angles. Note that the error ranges given in 
the plots of Fig.[3jb) are based on the uncertainty estimates of 
the fit routines. As additional uncertainties with regard to fit 
induced systematic errors are likely to exist, the true standard 
deviation for ^-factors of peak 2 may be significantly larger. 
Thus, it is not clear whether the different results for g-factors 
of peak 2 obtained for different 31 P concentrations are real or 
fit artifacts and therefore an unambiguous statement about the 
anisotropy of peak 2 is not possible. 

Previous pEDMR spectra on intrinsic c-Si have shown a 
second isotropic recombination signal with g 2.0023(6) 31 . 
This is in good agreement with peak 2 observed on c-Si sam- 
ples with [P] = 10 15 cm~ 3 . As with the previous study on 
intrinsic silicon, it is difficult to unambiguously assign peak 
2 to a particular interface defect type. Due to its strong in- 
homogeneity (the large line width) and its overall magnitude 
in comparison to the other pEDMR lines, it is possible that 
peak 2 is due to one or more randomly oriented anisotropic 
centers whose powder spectra would cause a macroscopically 
isotropic wide resonance line. The g-factor and the line width 
are close (yet not equal) to the g-factor and line width of sil- 
icon dangling bonds in amorphous SiC>2 (so called E' cen- 
ters^). Thus, it is possible that the observed pEDMR data 
is due to transitions involving E' centers in proximity of the 
c-Si:P/SiC>2 interface. Note that there are several distinct E' 
centers which differ by the backbonds of the Si atom on which 
the center is located (one or more oxygen atoms, hydrogen 
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atoms) or their relaxation state (the unrelaxed E' center has 
been called E' g ) 25 . Charge carrier trapping and recombination 
may work similarly for both the Pb centers with well defined 
orientation as well as the randomly oriented E'. Due to the 
localization length of the 31 P donor electron state 35 , it is also 
conceivable that transitions between the shallow donor states 
and the deep interface states are possible. Hence, peak 2 is as- 
signed in the following to an E' center with unspecified nature 
except for its proximity to the c-Si:P/Si02 interface. 



C. Identification of P to interface— defect transitions 

1. Experimental data 

The plots of Fig. [3] demonstrate the involvement of at least 
three qualitatively different paramagnetic electron states ( 31 P , 
Pb , E') in spin-dependent charge carrier recombination tran- 
sitions. This observation, however, does not prove that the 
observed recombination transitions take place between these 
different states. While the spin-dependency of transitions 
through localized states requires the existence of pairs of 
paramagnetic states 44 , there are examples of spin-dependent 
transitions which produce only a single resonance line in 
EDMR experiments 28,39 when transitions occur between iden- 
tical centers or when the paramagnetic pairs are strongly cou- 
pled. Therefore, the detection of multiple EDMR lines (as 
in Figs. |2] and [3J leaves the exact nature of these transitions 
elusive. 

One approach to an understanding of whether the observed 
^-factors belong to the same or different transitions is to ana- 
lyze the dynamics of the spin-dependent processes associated 
with these different resonances: Spin-selection rules usually 
discriminate permutation symmetries of spin s = 5 pair a 40,42 
which means that it is the mutual orientation of the two spins, 
not the individual spin state of one of the two pair partners, 
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Figure 4: Plots of the zero crossing times Tq of magnetic resonantly 
induced photocurrent transients of the 31 P low or high field resonance 
(Tq ) versus the zero crossing times of magnetic resonantly induced 
current transients of the Pb /E' center resonances (Tq ) measured for a 
variety of different samples and sample conditions (temperature, off- 
set current, illumination). The solid line is a linear function through 
the origin with slope 1 . The two plots made for the two concentra- 
tions show a strong correlation of the zero crossing times of all three 
resonances. 



which determines the transition rate. Therefore, the transient 
behavior of spin-dependent transition rates exhibits an identi- 
cal behavior after a spin resonant manipulation of either one of 
the two pair partners. If two EDMR detected resonances ex- 
hibit a different transient behavior after the same pulsed exci- 
tation, the spin-dependent transitions corresponding to these 
resonances must be different as well. However, if the transient 
behavior is identical, the two resonances may belong to iden- 
tical transitions. While different dynamics of spin-dependent 
transition rates for different ^-factors is proof that they belong 
to different processes, identical transition rates only indicate 
that the two resonances may be due to the the same process, 
as they may also be due to different processes which coinci- 
dentally have the same transient behavior. Therefore, the ob- 
servation of identical transients requires further testing of the 
correlation of the observed dynamical behavior under varying 
experimental conditions. 

Figure |2|b) demonstrates that the fast relaxing current 
quenching and the slowly relaxing current enhancing behav- 
ior that is discussed above occurs similarly for both the 31 P 
hyperfine resonances and the resonances associated with in- 
terface states. The dynamics (which means the decay times 
of the different exponential functions), but not necessarily 
the absolute magnitude of current transients from different 
centers involved in the same electronic process are identical. 
Thus a comparison of different transients yields information 
about the paramagnetic centers between which transitions take 
place. We compare the "zero-crossing time", To, defined as 
the time after the pulsed excitation of a spin-resonance in- 
duced current transient when the quenching and enhancement 
are identical. The comparison as shown in Fig. |2jb) clearly 
reveals identical To rs 18 jJ,s with an error of s» 1 jj,s. This 
suggests that the processes connected to these resonances are 
due to transitions involving both the 31 P donor states as well 
as the Pb interface state or the E' near-interface state in the 
way depicted in the sketch of Fig.QJb). 

In order to test whether To for the interface defects and 
the 31 P remained identical when the dynamics of the spin- 
dependent current signal is changed, the experiment presented 
in Fig.[2{b) was repeated under various combinations of tem- 
peratures (5 K, 8 K, 10 K, 13 K and 15 K), sample currents 
(10 jUA ~ 300 [lA (and therefore different electric fields), 
and sample surface orientations for the two different 31 P dop- 
ing concentrations mentioned above. Due to the variation 
in charge carrier concentration as well as Fermi- and quasi- 
Fermi energies caused by changing these experimental param- 
eters, the dynamics of the observed spin-dependent recombi- 
nation transitions, and therefore the dynamics of the observed 
current signals, changed significantly between measurements. 
For all data sets, To was determined for transients recorded at 
magnetic fields corresponding to the two 31 P hyperfine reso- 
nances as well as at the maximum of the overlapping interface 
defect signals. The results of this procedure are displayed in 
Fig.|4]in two correlation graphs, for samples with the two dif- 
fernet donor concentrations. Each graph displays a plot of Tq 
versus Tq with A and B corresponding to the 31 P low field and 
the 31 Phighfield resonances, respectively (represented by the 
circles) the 31 P low field and the interface state resonances, re- 
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spectively (represented by the squares) and the 31 P high field 
and the interface state resonances, respectively (represented 
by the triangles). 



2. Discussion 

The data in Fig.|4]reveals three observations: 

(i) As anticipated, the variation of experimental parameters 
clearly varies the dynamics of the observed signals. 

(ii) While similar combinations of experimental parameters 
were used for the two concentrations, the times To of the sam- 
ples with [P] = 10 15 cm~ 3 are generally slower than the To 
obtained for [P] = 10 16 cm~ 3 . This observation shall not be 
discussed quantitatively in the following, however, it is con- 
ceivable that as an increase of the 31 P doping concentration 
leads to a smaller main-pair distance (note the definition of 
main-pair distance given in section ITV A 21 for 31 P to inter- 
face defect pairs as well as 31 P to 31 P pairs. It should not lead 
to a change of the main-pair distance of interface defect pairs. 
The observed general decrease of the transition time with an 
increase of the 31 P concentration therefore suggests that the 
observed signals are predominantly due to transitions involv- 
ing 31 P and not transitions between two interface defects. We 
note that there may be other reasons for this decrease in transi- 
tion times with increasing doping density: one such possibil- 
ity is that different internal fields may lead to different energy 
separations between pair partners with different intra-pair dis- 
tance, with a corresponding change in transition probability. 

(iii) The entire set of data reveals a strong correlation of To 
between any combination of signals, whether it is between the 
two 31 P hyperfine signals or between the interface signals and 
either one of the two 31 P hyperfine signals. The correlation 
between the two 31 P hyperfine peaks is expected as the only 
difference between 31 P atoms contributing to these two peaks 
is the nuclear spin state, which has little influence on elec- 
tronic transitions times. The strong correlation between any 
two of the 31 P peaks and the interface states is again strong 
evidence that the spin-dependent transitions measured involve 
31 P -interface defect pairs in the way sketched in Fig. [Tib). 

It was not possible to obtain unambiguous correlation plots 
similar to those in Fig. [4] for the two different interface and 
near-interface defects discussed above. It is therefore not 
possible to verify whether the observed correlation between 
the 31 P and the interface defects applies to both interface cen- 
ters or only to the one which dominated under the observed 
conditions. However, we point out that the strong correlation 
seen in the plots of Fig. [4] is consistent with the assumption 
that the dynamics of both defects correlates with the 31 P dy- 
namics and thus, spin-dependent recombination through both 
interface defects in the manner sketched in Fig. [TJb) seems 
possible. 

The correlation data of Fig. |4] raises a question about 
the presence of spin-dependent transitions between inter- 
face defects. It is known from previous pEDMR studies of 
(111) oriented nominally intrinsic c-Si/SiC>2 interfaces at X- 
Band 29 - 30 - 31 , as well as from recent pEDMR studies at high 
magnetic fields (m 8.5T)2£, that spin-dependent transitions 



involving only interface defects but not 31 P are possible and 
that these transitions can be due to spin-dependent resonant 
tunnelling between two interface states in sufficient spatial 
and energetic proximity. We conclude from the data seen in 
Fig.|4]that, whilst such interface defect-only processes may or 
may not have been present at the investigated c-Si:P/SiC>2 in- 
terfaces, they did not dominate the observed spin-dependent 
rates, consistent with the very weak signal strength of previous 
pEDMR measurements on intrinsic c-Si/SiC>2 interfaces con- 
ducted at X-Band 31 . The interface defect signals are weak in 
spite of the presence of a significantly higher interface defect 
density compared to the areal density of 31 P close to the inter- 
face. This may be explained by considering the signal from 
interface defect pairs, which have almost identical g-factors 
and therefore stronger coupling, leading to lower pEDMR sig- 
nals than those seen from 31 P -defect pairs which have quite 
different g-factors and weaker, but still finite, coupling —, and 
therefore dominate the signal. 

D. Electrical detection of spin coherence 

1. Experimental data 

The ability to perform spin-Rabi nutation is a crucial pre- 
requisite for the coherence time measurements using echo 
pulse sequences as it reveals the resonantly induced Rabi fre- 
quency for a given set of experimental conditions (e.g. for the 
applied microwave power) . The Rabi frequency is needed 
to determine the pulse lengths required to obtain the correct 
nutation angles during the pulse sequences. We performed 
transient nutation style experiments for the electrical detec- 
tion of spin-Rabi nutation. For these measurements, the pho- 
tocurrent change AI was integrated between two appropriately 
chosen integration times t\ and ti after the microwave pulse so 
we obtain a charge 

Q= f' 2 M{t)dt (1) 
Jti 

that is proportional to the number of spin-dependent tran- 
sitions induced by the resonant excitation 42 . As Q is pro- 
portional to the projection of the resonantly prepared coher- 
ent spin state |v(t)) at the end of the excitation pulse of 
length t onto the singlet state \S) 10l4 ° , the measurement of 
<2(t) « \ (S\ v(t)) | 2 reveals the propagation of the singlet con- 
tent of |v(t)) during the pulse in a transient nutation-style 
experiment. Figure[51a) displays the measurement of Q(t,Bq) 
as function of the pulse length t and the magnetic field Bo for 
a c-Si sample with [P] = 10 16 cm~ 3 , a temperature T = 5 K, 
9 = 90°, and integration times t\ = 6 jus and tj_ = 16 /is. 

2. Discussion 

The data set displayed in Fig. [5] shows that Q(t) has an 
oscillating behavior around the same magnetic fields which 
produced local current response maxima in the data set dis- 
played in fig. |2] The oscillatory dependence on T is due to the 
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dephasing spin-Rabi nutations as demonstrated, for the mag- 
netic field Bo = 348 mT, by the data sets displayed in Fig.|5jb): 
The three plots show Q(l) for three different microwave pow- 
ers (Z?i field strengths). We anticipate 40 the spin-Rabi nuta- 
tion signal of an inhomogeneously broadened spin ensemble 
to follow the integral of a first kind Bessel function 

r rB,t i roo sin 2 (yBixVl+x 2 ) 
GO) / M2x)dx = - / ^ = '-dx 

(2) 

in which y w 2.8 MHz is the gyromagnetic ratio. A fit of 
the three data sets in Fig. [5jb) with Equ. [2] shows a good 
agreement and it reveals values for the nutation frequencies 
£lg = yB\ for the three different powers. A plot of Q.r versus 
the B\ fields obtained from the relative microwave attenuation 
used for the three measurements shows that the expected lin- 
earity of the spin-Rabi nutation is given. The data shown in 
Fig.|5]is thus shown to be due to the spin-Rabi nutation of 31 P 

x(ns) 




B i(mT) 

Figure 5: (color online) (a) Plot of the measured integrated charge 
<2(t,Bo) as defined by Equ. [T] as a function of the magnetic field 
Bq and the length T of an applied microwave pulse with frequency 
/ = 9.7475 GHz and power P = 8W. (b) The symbols represent a 
plot of the measured charge <2(t) for four different microwave pow- 
ers at Bq = 347 mT. Note that for the latter case, the pulse length was 
recorded up to T = 120 ns only as the measurement was limited by 
signal perturbation due to the pulse induced microwave current arti- 
facts. The blue lines represent fits of an integrated Bessel function 
to the experimental data. The Rabi-nutation frequencies obtained 
from these fits are displayed in the inset as a function of the applied 
B\ field. The fit of a linear function through the origin (black line) 
shows good agreement. 



donor and interface electron spins, confirming the previously 
reported observation made for c-Si(100):P/SiO2 interface with 
higher 31 P concentration. Note that with the data obtained 
from the transient nutation measurement, it is possible to de- 
termine the length of %- and ^-pulses as needed for the Ti 
measurements in the following. The latter will be necessary 
since the decay of the Rabi nutation as displayed in Fig. [5] 
is not a measure for spin-coherence. The agreement of the 
nutation data with the integrated Bessel function is indicative 
of coherent dephasing, not coherence decay, being the domi- 
nant source of the observed nutation decay. This assumption 
is confirmed by the electrically and pEPR detected echo data 
discussed below which shows that the real T2 spin coherence 
time of the 31 P donor electrons is distinctively longer than the 
decay of the nutation signal. 



E. Comparison of the coherence time T2 of P donor electrons 
at the c-Si:P/SiC>2 interface and in the c-Si:P bulk 

1. PEDMR and pEPR detected spin-echoes 

The data presented in sections fl V Al and II V CI showed that 
the spin-dependent current observed at c-Si:P/SiC»2 interfaces 
are due to transitions that take place between 31 P donor elec- 
tron states and c-Si:P/Si02 interface states in their proximity. 
A measurement of the 31 P donor electron spin coherence time 
(the T2 time) using pEDMR will therefore reveal only the T% 
of 31 P donor electrons in close proximity of the c-Si:P/Si02 
interface defects. In contrast, a 7*2 measurement using conven- 
tional radiation detected pEPR will reveal the coherence time 
of 31 P donors in the bulk as the contribution of near-surface 
31 P atoms in the pEPR signal will be negligible in comparison 
to the magnitude of the bulk 31 P signal. Thus, the compar- 
ison of T2 times measured by pEDMR and pEPR allows the 
influence of interface defects on the 31 P donor electron spin 
coherence times to be studied. 

The pEPR and pEDMR experiments require different ap- 
proaches to the measurement of T2 times: With pEPR, the T2 
times of paramagnetic centers can be determined most easily 
using a Hahn-echo pulse sequence. Hahn-echoes are a tem- 
porary rephasing of a spin-ensemble due to a pulse sequence 
consisting of an initial 7r/2-pulse which turns the ensemble 
polarization into the plane perpendicular to the Bo field and a 
subsequent 7£-pulse which initiates the phase reversal 4 ^. As 
illustrated in the inset sketch of Fig. [6jb), when the echo pulse 
sequence consists of the § — % pulses with a duration z be- 
tween the pulses, a Hahn-echo can be observed at a time T 
after the second pulse which is the time 2t after he first pulse. 
When a Hahn-echo is observed, T^-times can be measured by 
determining the decay of the Hahn-echo as a function of twice 
the pulse separation time 2t. Figure|6tb) shows a data set for 
a Hahn-echo transient obtained from a 31 P doped c-Si sample 
with [ 31 P] = 10 16 cm~ 3 at a temperature T = 10 K and with 
9 — 90°. One can clearly see a local maximum of the transient 
microwave signal at a time t = 2t. The data set was well fit 
by a Gaussian function which was used in order to determine 
the integrated intensity of the echo. 
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In contrast to pEPR measurements, pEDMR does not al- 
low direct observation of Hahn-spin echoes through real time 
transient measurements as the integrated sample current Q{i) 
always represents a projection of the spin-state at the end 
of the pulse sequence onto a singlet state. Therefore, the T2 
time measurement using pEDMR requires the utilization of a 
modified Hahn-echo pulse sequence that is illustrated in the 
inset of Fig. |6|a) 6 . In order to resolve the dynamics of the 
spin-ensemble during and after the ^ — ?r pulse sequence, a 
third pulse with length f is applied at a time z' after the be- 
gin of a conventional Hahn-echo pulse sequence. The third 
pulse projects the spin-ensemble at the time t' onto the in- 
direction which in turn determines the singlet content of the 
31 P -interface defect pair. The charge Q integrated following 
this pulse therefore represents the polarization of the spin- 
ensemble along the x- axis [y-axis] at the time t' (assuming 
the B\ field is directed along the y- axis [x- axis]). The mea- 
surement of the entire echo transient using the pEDMR detec- 
tion requires repetition of the echo sequence for various t'\ the 
third pulse (also called the detection pulse) is swept through 
the time range which covers the Hahn-echo maximum as well 
as the echo baseline. While this procedure makes the time 
needed for T2 measurements significantly longer, it allows the 
measurement of a T2 using pEDMR. The plot displayed in 
Fig. |6j a) shows an electrically detected spin-echo measured 
on the same sample and under identical conditions, recorded 
during the same experimental run as the measurements shown 
in Fig. |6jb). The data clearly shows the echo whose fit with 
a Gaussian function is displayed by the solid line. The com- 
parison of the two echo functions shows that the electrically 
detected echo is more narrow than the radiation detected indi- 
cating that it is due to a more heterogeneous spin ensemble. 

The use of two different measurement techniques (pEDMR, 
pEPR) raises the question of whether both methods probe the 
same observable, namely the T2 of 31 P impurity atoms. There 
has recently been a comparative study of T2 times confirming 
this identity using a pEDMR detected spin-dependent bulk 
process 37 (a spin-trap process of 31 P in c-Si that becomes rele- 
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Figure 6: Comparison of (a) an electrically detected spin-echo with 
T = 300 ns and (b) a conventional radiation detected Hahn-echo with 
% = 10 fls. The data sets were recorded on the same sample, under 
identical sample conditions (T= 10 K, / = 250 /lA, 6 = 90°) during 
the same experimental run. Both data sets were fit with Gaussian 
functions (solid line). Insets of the two plots show sketched timelines 
of the used pulse sequences. 



vant at high magnetic fields) which showed that both pEDMR 
and pEPR measured T2 times reveal an excellent agreement. 
Thus, in the following, systematic measurements comparing 
pEPR detected ^-times of 31 P bulk impurities and pEDMR 
detected ^-times of 31 P interface impurities are presented. 
These measurements are made for three reasons: (i) To ex- 
tend the previous observation of electrically detected Hahn- 
echoes 6 to c-Si(lll) surfaces, (ii) To measure the tempera- 
ture dependence of the 72-times. (iii) To obtain comparable 
measurements with both pESR and pEDMR under identical 
conditions. 



2. Measurement of spin-echo decays 

The spin-echo effects shown in Fig. [6] are imprints of co- 
herent spin motion on currents or radiation intensities. When 
a spin ensemble looses coherence during a -| — % sequence, 
the intensity of the spin-echo following this sequence decays. 
Quantifying the decay of Hahn-echoes is a direct measure of 
the coherence time 72 45 . Figure|7]displays a set of electrically 
detected Hahn echoes as well as a plot of the integrated echo 
intensities obtained from them, recorded on a sample with 
[P] = 10 16 cm~ 3 at magnetic fields in resonance with the high 
field 31 P peak at T = 10 K and = 90°, and with a sample 
current of / = 250 jiA. One can clearly see the gradual decay 
of the echo intensity with increasing pulse separation time t. 
The solid line of the plot in Fig.|7]displays a fit of the intensity 
data with a modified exponential function 

I(2T)=e 72 ^ (3) 

which contains a contribution due to a single exponential T2 
decay as well as the stretched exponential contribution due 
to the isotopical influence of the 29 Si which causes spin- 
diffusion with time constant Tg^2-. The plot in Fig. [7] ex- 
hibits a good agreement of the echo decay with the fit func- 
tion and thus, by using this method for both pEDMR as well 
as pEPR detected echo decay measurements, the T2 times of 
31 P impurities can be determined for the c-Si:P bulk and the 
c-Si:P/SiC>2 interface, respectively. Note that all pESR de- 
tected and most pEDMR detected echo measurements were 
conducted on the low field peak of the hyperfine split 31 P res- 
onance. In order to confirm that the spectral proximity of the 
Pb -peak to the low field 31 P resonance under the given ex- 
perimental conditions does not distort the electrical T2 mea- 
surements, control measurements were carried out on the high 
field 31 P peak for T = 5K. Within the given error margins, 
the results of these control measurements (Fig. [8] blue solid 
diamonds) are in agreement with the measurements obtained 
from the low field 31 P resonance. 



3. Temperature dependence of T2 times 

The measurement of T2 as described above was repeated on 
the same sample with pEPR at temperatures T — 5 K, 8 K, 10 
K, 1 1 K, 12 K, 13 K, 15 K and with pEDMR at temperatures 
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Figure 7: Plot of integrated intensities of the electrically detected 
echoes as a function of 2t recorded with pEDMR on a c-Si:P/SiC>2 
sample with [P] = 10 16 cm~ 3 . The solid line represents a fit with a 
modified multiexponential decay function (see text), the small inset 
plots represent plots of echo data sets. 



Figure 8: (color online) Plot of the normalized integrated echo inten- 
sity of pEPR and pEDMR detected Hahn echoes as a function of the 
logarithms of the pulse separation time T for various temperatures. 
The data sets were fit with the modified exponential function given 
by Equ.[3] All pEDMR and pEPR measurements were conducted on 
the same samples during the sample experimental run. 



of T = 5 K, 10 K, 13 K. The integrated echo intensities of 
these measurements are plotted in Fig.|8]along with the results 
of their fit with Eqn.[3] The comparison of the pEDMR and the 
pEPR data sets shows that while there is a strong temperature 
dependence of the echo decay for the pEPR data, the pEDMR 
decay is faster and, within the range T = 5 K to 13 K, nearly 
constant. In order to analyze these observations quantitatively, 
the fit results for T2 are plotted for both the pEDMR and the 
pEPR measurements in an Arhenius plot in Fig. [9] 

The coherence time of the bulk 31 P donors, T%, determined 
via pEPR, is well fit with a function of the form: 



1 =M 1 

— =AeV H 

T 2 T 



(4) 



where AE is an activation energy, and 1/Tq a constant relax- 
ation rate independent of temperature. The best fit occurs with 
AE = 9.1 ± 0.5 meV and l/T = 10 ± 1 kHz. 
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Figure 9: Plot of the logarithm of coherence decay rate T^ 1 obtained 
from the data shown in Fig.[8]as a function of the inverse temperature 
T . The solid lines are fits of the data. The pEDMR data is fit with 
a constant function. The pEPR is fit with a combination of a constant 
function and a temperature activated Boltzmann factor. 



As the pEDMR and the pEPR data displayed if Fig. [9] was 
collected on the same sample in the same experimental run, 
it is clear that the pEDMR measured T2 times of 31 P donor 
spins in close proximity to interface defects are significantly 
shorter than the pEPR measured T2 times and, within the er- 
ror, independent of the temperature. At the same time, the 
pEPR measured T2 of the bulk 31 P donor spins is not only 
longer than the pEDMR measured 31 P donor spins, it is also 
highly temperature dependent. In the temperature range from 
T = 5 K to 15 K, the coherence time changes by about a 
factor 50, exhibiting a good agreement with the previously 
demonstrated thermal activation 7 of 31 P in a c-Si environment 
with a natural abundance of the 29 Si isotope. The compari- 
son of T2 of 31 P in proximity to interface defects with those 
in the bulk therefore shows that interface defects significantly 



shorten the donor electron spin coherence time. Within the 
given temperature range, T2 appear to be pinned at as 1.3 jJ.s, 
a value which has been observed previously for electrically 
detected 31 P spins in c-Si samples with different surface orien- 
tations, donor concentrations and experimental conditions'^. 
The independence of the donor spin coherence time of near- 
interface defect 31 P atoms from experimental conditions, in- 
cluding temperature, suggests that in contrast to bulk donors, 
T2 of the near-surface donors is not determined by spin-spin 
or spin-lattice interactions but by transition directly related to 
the interface defects in their immediate proximity. 

We identify two possible origins for the drastic quenching 
of the 3 1 P donor electron spin T2 time in proximity of interface 
defects described above: 
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(i) The electronic transition between the 31 P donor state and 
the Pb center. This occurs when the donor electron falls into 
the doubly occupied interface ground state. The electronic 
transition leaves both the 31 P donor and the interface state dia- 
magnetic since there is no donor spin present after the transi- 
tion and the interface state is a doubly occupied singlet state. 

(ii) A fast T2 of the interface state which significantly 
quenches the Ti relaxation of 31 P donor spins in its proximity 
as suggested by de Sousai^. 

Note that the electronic transition is a limitation to both T\ 
as well as T% processes of the spin pairs as it destroys the pairs 
and when the electronic transition is the limiting factor for 
the T2 times, both, the electrically measured T2 and T\ times 
should be equal. In contrast, it is conceivable that the interface 
state induced enhancement of the 31 P Tj relaxation keeps the 
T\ either unchanged or significantly slower than the T2 time. 
As a consequence, if the electronic transition time is longer 
than the Pb induced T2 time, the electrically measured T\ time 
equals the electronic transition time and it should therefore be 
significantly longer than the electrically measured T2 times. 
Hence, a study of T\ times of the near-interface defect 31 P 
donor electron spins, as presented in the following, has been 
conducted. 
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Figure 10: Demonstration of an electrically detected spin inversion 
recovery experiment. The inset is a sketch of the inversion recovery 
pulse sequence which consists of the Hahn-echo sequence that is 
preceded by an inversion (n-) pulse at a time t" before the Hahn- 
echo sequence begins. Similar to the electrically detected Hahn-echo 
shown in Fig. [6] a projection pulse is shifted through the sequence 
during different repetitions of the experiment. The main plot displays 
an echo recorded by plotting Q as a function of the difference t' — 2t 
between the projection pulse begin t' and the echo maximum at 2t. 



2. Inversion recovery of near-interface defect P donor electrons 



Comparison of T2 and the longitudinal relaxation time T\ of 
31 P donor electrons at the c-Si:P/Si02 interface 



1. Electrical detection of spin inversion 

In order to probe T\ of near-interface 31 P , electrically de- 
tected inversion recovery experiments were carried out for a 
temperature of T — 5K where pESR measurements of bulk 
31 P donor spins reveal large differences between T\ and T2 
times 7 . The idea behind the inversion recovery experiment 45 
is to invert a given spin polarization towards a negative non- 
equilibrium polarization and to then observe transiently the 
gradual return toward the equilibrium due to T\ processes. 
Thus, the experiment consists of an initial inversion pulse, a 
defined delay time t" during which the spin ensemble under- 
goes T\ relaxation and a subsequent polarization measurement 
which is typically conducted by utilization of a Hahn-echo 
sequence with fixed pulse separation time t. Figure [10] dis- 
plays a sketch of the pulse sequence used for the electrical 
inversion recovery measurements presented in the following. 
Similar to the Hahn-echo decay measurements, the Hahn- 
echo is measured with pEDMR by repetition of the experi- 
ment whilst applying projection pulses which are gradually 
swept through the echo sequence. The data in Fig. [TOl shows 
an inverted Hahn-echo recorded with a short t = 252/xs and 
an equally short delay time t" = 52ns producing strong inver- 
sion. Instead of the positive spin echo as detected for a simple 
Hahn-echo sequence, the sign of the echo is negative. To the 
knowledge of the authors, the data shown in Fig. [10] is the 
first demonstration of an electrically detected spin inversion 
experiment. 



In order to determine T\, the experiment shown in Fig. [TOl 
was repeated for seven values of t" in the range of 52ns < 
t" < 20/Xs. The results of these experiments are displayed in 
Fig- HU They show that the polarization inversion that exists 
directly after the inversion pulse (t" is very small) exponen- 
tially approaches the steady state polarization with increasing 
t". The integrated echo amplitudes were obtained from a fit 
of the echo data with Gaussian functions. They dependence 
on t" show an excellent agreement with an exponential decay 
function with a negative offset 



M(t") =M(0) 



\-2e~ 



(5) 



. The time constant T\ = 4.0(5 )jj.s obtained from this fit is 
more than six orders of magnitude shorter than the previously 
investigated bulk T\ times 7 which shows that the proximity of 
Pb centers leads to dramatically reduced T\ times. 



3. Discussion 

The measurement of T\ of 31 P donor electrons in proximity 
of interface defects revealed a T\ time which was quenched 
by orders of magnitude compared to bulk 31 P : However, in 
contrast to bulk pESR measurements, T\ is only about a factor 
3 longer than the electrically detected T2 time at the measured 
temperature of T = 5K. Following the argumentation given 
in section II V E 41 we learn that since T2 < T\ , T\ will be de- 
termined by the electronic transition but not T2 as it is faster 
due to the influence of the nearby interface defect. It shall be 
emphasized though that the errors of the T\ and T2 measure- 
ments have been deduced from the fits of the Hahn-echo and 
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Figure 11: (color online) Plot of the integrated echo intensity de- 
tected with the pulse sequence shown in Fig.[l0]as a function of the 
inversion recovery time t". The solid line represents a fit of the data 
with a single exponential function. Note that the inverted echo for 
small t" changes into a non-inverted echo with equal magnitude for 
large t". The inset plots on the left show the raw data of the various 
echo measurements as well as fits with Gaussian functions which are 
the basis for the integrated echo intensities. 

inversion recovery data with eqs.[3]and|5] respectively which 
leaves a significant uncertainty regarding the magnitude of the 
measured T\ and T2 times. Given the proximity of the two val- 
ues within the calculated errors and the remaining uncertainty 
about the given errors, it can not be absolutely excluded that 
7] w T2. The latter would be indicative that the electronic tran- 
sition not only determines T\ but also ?2. Note that in spite of 
this uncertainty, is is clear that the electronic transition can not 
be slower than the measured Ti-time and it can not be faster 
than the measured ^-time. Hence it is unambiguously clear 
that the 31 P -interface defect process takes place on a time 
scale in the lower microsecond range (Ifls — 4/xs). 



V. CONCLUSIONS 

The realization that T\ and T2 times of 31 P donor electron 
spins near-c-Si:P/SiC>2 interface defects are drastically short- 
ened to a similar value is proof that the recombination of the 
31 P donor electron into a charged interface ground state is the 
lifetime limiting process of the 31 P qubit. This conclusion is 
also supported by the independence of the transition time on 
the temperature as this suggests that the observed spin relax- 
ation times are determined by the 31 P -Pb pair geometry and 
not phonon density. This conclusion is also consistent with the 
assumption that only pairs with a certain, narrowly distributed 
"main" distance contribute to the observed pEDMR signals. 



Our results have implications for proposed mechanisms to 
readout donor spins using nearby probe spins, particularly in- 
terface defects. Without the ability to vary the coupling be- 
tween donor and probe spins, the electronic transition times 
will be fixed, limiting T2 and also limiting the readout time. To 
overcome this problem, control of the coupling between target 
and probe spins will be required. This may be achievable us- 
ing the Stark effect as electric fields will have a stronger effect 
on the localization of the slightly delocalized 31 P donor wave 
function in comparison to the extraordinary strongly localized 
deep interface defects. Whether these different electric field 
sensitivities have a significant impact on the exchange cou- 
pling between the 3 1 P donor interface states in close proximity 
at sufficiently low fields remains to be proven. 

The data presented and discussed above strongly supports 
the model for spin-dependent recombination via 31 P donors 
and interface defects that was presented by Stegner et alJ£. 
Whilst this explanation for the observed EDMR signal has 
become commonly accepted, and is further supported by the 
results presented here, the possibility remains that other pro- 
cesses also contribute to the resonant changes in current. 
Given the now significant evidence supporting the 31 P-Pb 
model 6 ' 1Q i 27 , we agree with the conventional understanding 
and conclude here that we are indeed observing spin depen- 
dent transitions between 31 P donors located close to Pb defects 
at the Si(l 1 1)-Si02 interface. 

Finally, we note that the effect of spin-dependent transi- 
tions in reducing coherence times may not be restricted to 
the mechanisms presented here, but may also be important for 
other mechanisms, such as spin dependent scattering of con- 
duction electrons by donor spins in two dimensional electron 
gasses 4VW>r*>. 



VI. SUMMARY 

In summary, we have reported on the investigation of 
spin-dependent processes at the c-Si:P/Si02 interface using 
pEDMR, and shown that spin-dependent 31 P to interface de- 
fect recombination takes place at the c-Si (111) surface in a 
similar way as previously described processes for c-Si(100) 
surfaces. The imprints of spin-dependent recombination on 
interface currents reveal EPR resonances of the hyperfine split 
31 P resonance and Pb defects and at least one other defect, 
assigned here to the unrelaxed E' defect in the Si02 matrix, 
previously unobserved in 31 P doped samples. The correla- 
tion measurements of the dynamics of these pEDMR detected 
signals strongly support the model that recombination tran- 
sitions between 31 P and the interface defect states occur. By 
electrical detection of spin-echoes, a measurement of both the 
coherence time T2 as well as the longitudinal relaxation time 
T\ of 31 P donor spins in proximity to interface defects was 
possible. These measurements revealed that T2 ~ 1 .3 jtis inde- 
pendent of the applied temperature range for 5K <T< 13K. 
The longitudinal relaxation time T\ w 4/is at T = 5K. The ob- 
servations contrast the pESR measured T\ and T2 of bulk 31 P 
which are significantly longer and strongly temperature de- 
pendent. Our observation show that the proximity of interface 
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defects to 3 1 P donors introduces spin-dependent recombina- 
tion and that these electronics transitions significantly limit 
both, the phase coherence time T2 as well as the longitudinal 
relaxation time T\ of donors near interface defects. The impli- 
cations of these findings for possible applications of the 31 P 



-interface defect transition as spin 31 P readout for proposed 
potential spin electronics or quantum information applications 
have been discussed. 



* Electronic address: dane.mccamey@physics.utah.edu 
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